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Rare Kaon Decays 

Giuseppina AnzivincEl 
Department of Physics, University of Perugia, Perugia, Italy 

Recent results on charged rare kaon decays from the NA48/2 CoUaboration and on neutral kaon 
decays from the KLOE Collaboration will be reviewed. 
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I. INTRODUCTION 

The results presented in this paper are mainly pre- 
cision tests of Chiral Perturbation Theory (ChPT), 
with focus on radiative non-leptonic kaon decays from 
the NA48/2 and KLOE Collaborations. Semileptonic 
decays {Vus related), leptonic and LFV processes and 
future rare kaon decays experiments are covered by 
other talks in these proceedings. 



II. NA48/2 - CHARGED KAON DECAYS 

The NA48/2 experiment at CERN has taken data 
in 2003 and 2004 with the main purpose of searching 
for direct CP violation in the decays of charged kaons 
into three pions. However, the very high statistics 
collected, the world's largest, has allowed to study 
many other decays and to measure their Branching 
Ratios (BR). In this paper I will fucus on radiative 
nonleptonic decays, that provide a crucial test of the 
ChPT fill and can give information on the structure 
of hadronic interactions at low energy. 



A. The NA48/2 experiment 

Two simultaneous focused kaon beams of oppo- 
site charge, with a central momentum of 60 GeV/c 
and a momentum band of ±3.8% are produced by 
a 400 GeV proton beam impinging on a 40 cm Be 
target. The decay volume is a 114 m long vacuum 
tank; the final states are reconstructed by a mag- 
netic spectrometer and a liquid krypton calorime- 
ter (LKr). Charged particles are measured by the 
magnetic spectrometer, consisting of four drift cham- 
bers and a dipole magnet located between the sec- 
ond and the third chamber; the momentum resolu- 
tion is a{p)/p = 1.02% © 0.044%p (p in GeV/c). 
The magnetic spectrometer is followed by a scintil- 
lator hodoscope consisting of two planes segmented 
into horizontal and vertical strips and arranged in 
four quadrants (charged hodoscope). The electromag- 
netic Liquid Krypton calorimeter is an almost homo- 
geneous ionization chamber with an active volume of 



~ 10 nn? and a 27 Xq thickness; the energy resolution 
is a{E)/E = 0.032/\/;B©0.09/i;©0.0042 (E in GeV). 
The space resolution for single electromagnetic shower 
can be parametrized as (Jx — Oy — 0.42/-\/E©0.06 cm 
(E in GeV). At a depth of ^ 9.5Xo inside the active 
volume of the calorimeter, a hodoscope consisting of 
a plane of scintillating fibres is installed (neutral ho- 
doscope); the signals from the four quadrants are used 
to give a fast trigger. 

A more detailed description of the detector can be 
found in 0. 



B. The K- 



^e decay 



The — > TT^e+e^ decay is a FCNC process in- 
duced at one-loop level in the Standard Model and 
highly suppressed by the GIM mechanism. The dy- 
namics of the decay is completely specified by the 



invariant function W{z), where 



: {M^JMkY is 

a kinematic variable. Several models have been de- 
veloped predicting the form factors that character- 
ize the decay rate and the dilepton invariant mass 
distribution. In the present analysis the following 
parametrizations of the form factors are considered: 



• Linear: W{z) = GfM|/o(1 
normalization and slope {5z). 



5z). with free 



: W{z) = 

GFM%{a+ + b+z) + W'^iz) with free param- 



• Next-to-Leading Order ChPT 



*on behalf of the NA48/2 and KLOE Collaborations 



eters (a_|_,6_|_) and an explicitly calculated pion 
loop term W'^iz). 

• A ChPT model developed by a group from 
Dubna [1] involving meson form factors W{z) = 
W{Ma, Mp, z), with meson masses {Ma,Mp) 
treated as free parameters. 

The aim of the analysis is to extract the form fac- 
tor parameters in the framework of each of the above 
models and to measure the corresponding BR's in 
the full kinematic range and, in addition, a model- 
independent BR in the visible kinematic range {z > 
0.08). The K"^ 7r^e+e~ rate is measured rela- 
tively to ^^^D (with TT^ e"'"e~7). Since 
the two decays contain the same charged particles in 
the final state, common selection criteria have been 
used, resulting in cancelation of particle ID inefficien- 
cies at first order. At the end of the selection, based 
on the reconstruction of three-track events, a total 
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of 7146 7r±e+e- candidates with 0.6% back- 

ground are found in the signal region. The recon- 
structed (7r^e"'"e~) invariant mass spectrum is shown 
in Fig. m 




M(ji±e*e ), GeV/c^ 



distinguish between the models considered. The BR 
in the full kinematic range, which includes an uncer- 
tainty due to extrapolation into the inaccessible region 
z < 0.08, is 

BR = (3.08±0.04,t ±0.04,3,^.±0.08ea;t.±0.07„.)xlO"^ 
BR{K^ n^e+e-) = (3.08 ± 0.12) x 10"^ 

This result is in fair agreement with previous mea- 
surements. In particular, comparison to the most pre- 
cise BNL E865 result [5], dismissing correlated un- 
certainties due to external BR's and model depen- 
dence and using the same external input, reveals a 
1.4 CT difference. The measurement of 6 is in agree- 
ment with the previous measurements based on K'^ 
TT+e+e" d, Q and 7r*^+/z~ 01 samples, and 

further confirms the contradiction of the data to me- 
son dominance models The measured fo, a+ and 
6+ are in agreement with the only previous measure- 
ment 5]. The measured parameters Ma and Mp are 
a few % away from the nominal masses of the reso- 
nances 0. 



FIG. 1: Reconstructed spectrum of {tv e^e ) invariant 
mass: data (dots) and MC simulation (filled area). 

The computed values of dTT^^f,/dz vs z are shown in 
Fig. [5] with the results of the fits to the three consid- 
ered models; the measured parameters and the corre- 
sponding BR's are presented in Table HI 



C. The 7r*77 decay 
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FIG. 2: The computed dF-^ee/dz and the results of the fits 
according to the three considered models. 

Fits to all the three models are of reasonable qual- 
ity, however the linear form-factor model leads to the 
best (see Fig. [5]). The data sample is insufficient to 



The contributions of the chiral lagrangian to this 
decay [l3] appear at Oip^), where only the A/ = 
1/2 invariant amplitudes A{z) and C(2:), with z = 
M^^/M^±, contribute. The decay rate and the spec- 
trum strongly depend on a single parameter c pre- 
dicted to be positive and 0(1). The invariant M^^ 
distribution is favored above 2m^+ and exibits a cusp 
at 2m^+ threshold. Calculations at 0{p^) [ll| show 
that unitarity correction effects can increase the BR 
by 30 — 40%, while vector meson exchange contribu- 
tions would be negligible. The — > 7r^77 decay 
rate is measured relatively to the — > tt^tt*^ nor- 
malization channel; they have identical particle com- 
position in the final state and only the cut on the 77 
invariant mass differs for the two channels. About 
40% of the total NA48/2 data sample have been an- 
alyzed and 1164 candidates have been found with an 
estimated background of 3.3 % which has to be com- 
pared with the only previous experiment that 
collected 31 events. This decay and the normaliza- 
tion channel were collected through the neutral trigger 
chain, intended for the collection of — > Tr^Tr^Tr" 
and therefore suffered from a very low trigger effi- 
ciency (ss 50%). Elaborate studies were performed 
in order measure trigger efficiencies and correct for 
them. The reconstructed invariant 77 mass spectrum 
in the accessible kinematic region is shown in Fig. [31 
also shown is the MC expectation assuming the ChPT 
prediction of 0{p^) [iT|. The model dependent BR of 
— > 7r*77 has been measured assuming the va- 
hdity of the 0{p^) ChPT and c = 2. The following 
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TABLE I: Results of fits to tlie tliree models and the Model-Independent BR(z > 0.08). 



5 = 

/o = 
BRi X 10^ = 


2.35 ± O.lSstat. 
0.532 ± 0.0f2stat. 
3.02 ± 0.04,tat. 


± 0.09syst. 

± 0.008,y,t. 
± 0.04syst. 


± O.OOcxt. 
± 0.007cxt. 

± 0.08ext. 


= 2.35 ± 0.18 
= 0.532 ± 0.016 
= 3.02 ± 0.10 


a+ =- 
h+ 

BRa X 10^ = 


-0.579 ± 0.012etat. 

-0.798 ± 0.053stat. 

3.11 ± 0.04stat. 


± 0.008syst. 
± 0.037syst. 
± 0.04syst. 


± 0.007cxt. 
± 0.017oxt. 
± 0.08oxt. 


= -0.579 ± 0.016 
= -0.798 ± 0.067 
= 3.11 ± 0.10 


Mp/GeV = 
BRs X lO'^ = 


0.965 ± 0.028Btat. 

0.711 ± O.OlOstat. 

3.15 ± 0.04stat. 


± 0.018syst. 
± 0.007sy.t. 
± 0.04syst. 


± 0.002ext. 
± 0.002ext. 

± 0.08cxt. 


= 0.965 ±0.033 
= 0.711 ±0.013 
= 3.15 ± 0.10 


BRmi X lO'' = 


2.26 ± 0.03stat. 


± O.OSsyst. 


± 0.06cxt. 


= 2.26 ± 0.08 




0.2 0.22 0.24 0.26 0.28 0.3 0.32 0.34 0.36 

m.„ [GeV/c^] 



FIG. 3: Reconstructed invariant 77 mass [GeV/c^]. 

preliminary result has been obtained: 

BR{K^ ^ 7r±77) = (1.07±0.04,t,t.±0.08,j,,t.) x IQ-^ 

A model independent BR measurement is in prepa- 
ration, together with the extraction of c from a com- 
bined fit of the spectrum shape and the decay 
rate. 



GeV/c^ was applied. Using the full NA48/2 data sam- 



ple, 120 



TT e ' e 



'7 (6.1 % background, esti- 



mated by MC) candidates are found in the accessible 
kinematic region. 

Fig- HI shows the reconstructed spectrum of e"'"e~7, 
with the MC expectation for the background contri- 
butions. To determine the BR in a model independent 
way, a partial branching fraction has been computed 
for each 5 MeV/c^ wide TOe+e-7 interval i: 



TT^e'^e 



with N'!^"^^ and N^^^ the number of observed signal 
and estimated background events, A^*^^ the accep- 
tance in bin i. The overall trigger efficiency is e and 
$x is the total kaon flux. By summing over the bins 
above m^-^f,-^ ~ 0.26 GeV/c^, the model independent 
BR is obtained: 



7) = 



7ree7 



N' 



bkg 



A. 



7ree7 



K 



BR{K~- 



■ IT 



^e~^e--i){Mee^ > 0.26GeV/c^) = 
= (1.19±0.12, ' " ' 



2stat.±0msyst.) X 10-« 



The parameter c has also been measured assuming 
the vahdity oi 0{p^) The final results of the 

analysis have been recently published p^ . 



D. The if^ TT^e+e 7 decay 

The 7r*e+e~7 decay, observed for the first 

time by NA48/2, is similar to the tt^jj with 

one photon internally converting into a e'^e~ pair. Af- 
ter the basic selection criteria, the data were domi- 
nated by '^^'^D decay, where the 7r° under- 
went a Dalitz decay tt^ — > e^e^7. A number of ad- 
ditional selection criteria had to be applied in order 
to effectively suppress the remaining background due 
to decays and accidental activity. The signal re- 
gion was defined by requiring 480 < TO^ie+e-^ < 505 
[AleV/c^]. Since the ChPT predicts only small sig- 
nal rate and the background increases for low values 
of me+e-7, the additional requirement m^+^-^ > 0.26 




Invariant eVy mass [GeV/c^] 

FIG. 4: Reconstructed invariant e+e~7 mass [GeV/c^]. 
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III. KLOE - NEUTRAL KAON DECAYS 

The results presented in this section are based on 
data collected with the KLOE detector at DA$NE, 
the Frascati e+e~ collider operated at a center of mass 
energy of 1020 MeV, the mass of the (/)-meson, that is 
produced almost at rest and decays, in 34% of cases, in 
K^K^ pairs. The two kaons are always a KsKl pure 
pair, so detection of a Kl guarentees the presence of 
a Ks- This procedure, called tagging, provides a pure 
Ks beam. 



A. The KLOE detector 

The KLOE detector consists of a large cylindrical 
drift chamber 15] of 4 m diameter and 3.3 m length 
operated with a low Z and density gas, surrounded by 
a lead-scintillating fiber calorimeter (EMC)[lB|- The 
chamber provides tracking, measuring momenta with 
a resolution of S{p±/p±) of 0.4% at large angle and 
a vertex reconstruction resolution of ~ 3.3 mm. A 
superconducting coil around EMC provides a 0.52 T 
magnetic field. The EMC is ^ 15 Xq thick and covers 
98% of the solid angle. Energy and time resolutions 
are a{E)/E = 5.7%/Ve and a{t) = 57 ps/VE ® 
100 ps, respectively. 



B. The Ks 77 decay 

The Ks 77 decay is an important test of ChPT. 
The decay amplitude can be calculated unambigu- 
ously at the leading 0{p'^) order of the perturbative 
expansion with an uncertainty of only few percent [l7j , 
giving BR {Ks -> 77) = 2.1 x 10"*^. 

In the data sample analyzed, corresponding to an 
integrated luminosity of / Cdt - 1.9 fb'^, - 700x10^ 
Ks decays are identified by iCi-tagging. Because of 
the tagging, no background from Kl — > 77 decay is 
expected, the main background being Ks tt'^tt^, 
with the two photons undetected, due to geometrical 
acceptance or non reconstruction in the calorimeter. 
One of the most effective variables against this kind 
of background is the arrival time of the photons: a 
prompt photon is defined as a neutral cluster in the 
electromagnetic calorimeter satisfying the condition: 
|T — R/c\ < min{5 at, 2 ns), where T is the time of 
flight, R is the cluster position with respect to the de- 
tector origin of the coordinates and at the total time 
resolution. A signal-enriched sample is defined by re- 
quiring two and not more two prompt photons in the 
event. To maximize Ks —> tt'^tt'^ rejection, only clus- 
ters with E > 7 MeV and cos6 < 0.93 are considered; 
the residual background from Ks decays other than 
27r° is suppressed by vetoing events with photons ab- 
sorbed in a small angle calorimeter. The event count- 



ing is performed using a kinematic fit imposing seven 
constraints: energy and momentum conservation, the 
kaon mass and the two photon velocities. 

Two other variables with a powerful discrimina- 
tion against the background are the two photon in- 
variant mass, M-y-y, and the opening angle between 
the two photons in the Ks center of mass system, 
0*^. To obtain the number of Ks 77 events, a 
2 dimensional binned maximum-likelihood of the fi- 
nal sample distribution in the Mj^ and cosO*^ vari- 
ables is performed using the MC generated signal and 
background shapes, taking into account data and MC 
statistics. The result of the fit is: iV(77) = 711 ± 35, 
with x^/dof = 854/826, 24.3% CL. Projections of 
data and fit are shown in Fig. [5l The signal cos9*^ 
is peaked at cos9=-l, while My^ distribution is gaus- 
sian at the Ks mass. The background is less peaked 
at cos9=-l and shows a broader distribution in mass, 
populating low mass values. The BR is obtained 




FIG. 5; Distribution of cosO^j and Mj^ for the final sam- 



ple. 



from N{Ks — > 77) using as normalization channel the 
Ks — > tt^tt" events recorded in the same data sample 
by counting the Ks tagged events with four prompt 
photons. The number of Ks tt'^tt" events after cor- 
recting for selection efficiency is (190.5 ± 0.5) x 10^. 
Using the latest PDG 9] value for BR{Ks 7r°7r°) = 
(30.69 ± 0.05)%, the BR(ii's ^ 77) is obtained 

BR{Ks ^ 77) = (2.26 ± Q.l2,tat. ± 0.06,^,,*.) x 10-^ 

This result together with other existing measurements 
of BR{Ks 77) as weU as the ©(p"*) ChPT theoret- 
ical prediction are shown in Fig. [SI There is a 3a 
discrepancy between the KLOE result and the NA48 
latest measurement [18|: the NA48 measurement im- 
plies the existence of a sizable 0{p^) counterterm in 
ChPT, while the KLOE result makes this contribution 
practically negligible. 
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FIG. 6: Comparison of mi{Ks 
ChPT predictions. 



77) measurements and 



A signal box to select Ks e+e^ events can be 
conveniently defined in the — M^e plane. Other 
independent requirements, studied by MC and tuned 
in the sidebands, have been applied in order to reduce 
the background contamination, before applying the 
— Mee selection. The signal box is chosen with 
an optimization procedure based only on MC. The 
cut for the signal box definition has been chosen to 
remove all MC background events: < 70. The 
cut on Afee invariant mass is set at 477 < M^e < 
510 MeV/c^, which rules out all signal events with a 
radiated photon with energy greater than 20 M eV. 
The signal box selection on data gives Nobs = 0. The 
upper limit on BR(A'5 — s- e^e") is evaluated using a 
bayesian approach and normalizing to the number of 
Ks ■K'^TT~ events. The result is: 



The Ks — > e^e decay 



BR{Ks e+e-) < 9.3 x 10"^ (90% CL) 



The decay Ks — > e~*'e~, like Kl e'^e~ and Kl 
is a FCNC process suppressed in the SM and 
dominated by two photon intermediate state. Pre- 
cise theoretical predictions, based on ChPT at ©(p^), 
evaluate the ratio T{Ks e~^e~)/r{Ks 77) = 
8 X 10~^ with 10 % uncertainty [l^. Using the present 
average for BR{Ks 77) = (2.71 ± 0.06) x 10~^ 
0, the SM prediction is BR {Ks e+e") - lO^^^^. 
The best experimental limit has been obtained by the 
CPLEAR Collaboration BR {Ks e+e") < 

1.4 X 10"'' at 90 % CL. Ks decays are identified by 
-ftr^-tagging and selected by requiring two tracks of 
opposite charge with a good vertex. One of the main 
backgrounds is the Ks 7r"'"7r~ decay, with two pions 
misidentified; for these decays the Mee invariant mass 
is peaked at low values, so a cut at Mee > 420 MeV/c 
is effective in rejecting this type of background. An- 
other important background is the (f> — > 7r"'"7r~7r'^ de- 
cay with one prompt photon that simulates a K^ in- 
teraction in the EMC and the other photon is not 
detected. After preselection, the data sample is re- 
duced to ~ 10^ events. In order to improve the sepa- 
ration between signal and background, a x^-like vari- 
able is defined, using information from the clusters 
associated to the candidate electron tracks. Using 
the MC signal events, likelihood variables are built 
based on: the sum and the difference 6t of the two 
tracks, where 5t — td — L/ [3c is evaluated in the elec- 
tron hypothesis; the ratio E/p between the cluster 
energy and the track momentum, for both charges; 
the cluster position relative to the extrapolation of 
the track, for both charges. Fig. [7] shows the corre- 
lation between x^ and Mee for MC signal and back- 
ground. Two sidebands are defined to check the con- 
sistency of MC with data and normalization: region 
1 (Mee < 460 MeU/c^), dominated by Ks tt+tt" 
events, and region 3 (Mee > 530 MeV/c^), mostly 
populated by </) — *■ 7r"'"7r~7r'^. 



M 300 p 
250 - 
200 - 



100 - 
50 

" 450 500 550 600 650 700 

M^^ (MeV) 

FIG. 7: Distribution of as a function of the invariant 
mass Mee: MC signal (red), background events (black). 

This measurement improves by a factor ~ 15 the 
CPLEAR [2(?| result and includes for the first time 
radiative corrections. 

D. The Kl — > 7r*e^J/(7) decay 

The study of radiative Kl decays offers the possi- 
bility to obtain information on the kaon structure and 
to test ChPT. Two different components contribute to 
the photon emission: the inner bremsstrahlung (IB) 
and the direct emission (DE). In the K^y.-y decay the IB 
component is much larger than the DE and the inter- 
ference terms. It is customary to apply standard cuts 
to allow comparison among results: E* > 30 MeV and 
0* > 20 {E* and 9* energy and photon angle w.r.t. 
the Icpton in the kaon rest frame, respectively). The 
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ratio R is defined: 



R = 



T{K%^-E*^ > 30 MeV, 9* > 20) 



Theoretical prediction for R range between 0.95 
xlO-2 and 0.97 x lO'^ [p. The decay ampHtude 
can be written as 



dT dV 



IB 



dEl 



dEZ 



{x)f{E;) 



where the second term is the SD contribution. All the 
information on the SD term is contained in the effec- 
tive strength {X). Candidate Kl events are tagged 
by the presence of a Ks — > tt+tt" decay. To count 
K^^^ signal events, a two-dimensional fit in the vari- 
ables E*, 9* is performed; this allows to measure both 
R and (X) . The results of the fit are shown in Fig. [H 
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FIG. 9; Comparison of K and (X) results (see text). 

Taking into account all systematics, the measure- 
ment of R and (X) yields 



i?^ (924±23,t,t±16,j,,t) X 10- 



(X) = -2.3±1.3,tat±1.4 



FIG. 8: Fit and fit residuals in bins of E% 



syst 



The measured value of {X) is in agreement with 
0{p^) evaluation in [2l|. The presence of DE contribu- 
tion reduces the value of R of about 1 %. Fig. [9] shows 
the comparison of this results with previous measure- 
ments. The present accuracy is not sufficient to solve 
the experimental discrepancy between NA48 [l^l and 
KTeV m. 
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